Topological kinematic constraints: quantum dislocations and the glide principle. 
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Topological defects play an important role in physics of elastic media and liquid crystals. Their 
kinematics is determined by constraints of topological origin. An example is the glide motion of 
dislocations which has been extensively studied by metallurgists. In a recent theoretical study 
dealing with quantum dualities associated with the quantum melting of solids it was argued that 
these kinematic constraints play a central role in defining the quantum field theories of relevance to 
the description of quantum liquid crystalline states of the nematic type. This forms the motivation 
to analyze more thoroughly the climb constraints underlying the glide motions. In the setting 
of continuum field theory the climb constraint is equivalent to the condition that the density of 
constituent particles is vanishing and we derive a mathematical definition of this constraint which 
has a universal status. This makes possible to study the kinematics of dislocations in arbitrary 
space-time dimensions and as an example we analyze the restricted climb associated with edge 
dislocations in 3+ ID. Very generally, it can be shown that the climb constraint is equivalent to 
the condition that dislocations do not communicate with compressional stresses at long distances. 
However, the formalism makes possible to address the full non-linear theory of relevance to short 
distance behaviors where violations of the constraint become possible. 

PACS numbers: 61.72.Lk, ll.30.-j 



I. INTRODUCTION 



Topological defects in crystalsiaS* 3 ^* 5 - exhibit a rich va- 
riety of kinematic properties 6 * 7 * 8 -. Defects such as intersti- 
tials/ vacancies, dislocations, and disclination o 6 : 7 : 8 : 9 per- 
turb the ideal perfectly periodic crystal. To date, numer- 
ous works extensively studied the dynamics and coarsen- 
ing of these defects, e .gii2iiiii2iiiiiii^. In normal solids, 
dislocations are present at low concentrations and their 
peculiar 'glide motions' are an important factor in de- 
termining the plastic properties of the medium. It has 
long been recognized that the energy-entropy balance of 
topological defects is responsible (via deconfincmcnt) for 
melting transitiona 16 : 17 i 18 i 19 i 20 i 21 i 22,23 . In glasses, an ex- 
tensive configurational entropy of these defects 2 ^ (from 
which ensuing restrictive slow dynamics might follow^) 
may be sparked. Investigations of Frank Kasper phases 26 
of defect lines and of related systems 27 : 28 : 29 ! 30 : 31 : 32 : 33 
have flourished. In some solids, such as shape memory 
alloys 3 -^, reducing the appearance of defects is a matter 
of pertinence. The study of topological defects in liq- 
uid crystals**^, compounded by the very crisp imaging 
techniques of these systems, has evolved into a broad 
field. Topological defects often display sharp dynamical 
imprints. Perhaps the best studied such dynamical ef- 
fect is the glide principle which forms the focus of the 
current article. Throughout the years, much work has 
been carried out on the classical glide motion of disloca- 
tions. Many of these works entail detailed studies of the 
Peach-Kohler forces between dislocations and Peierls po- 
tentials in systems ranging from classical solids to vortex 



arrays to liquid crystals, e.g. 36 : 37 : 38 : 39 : 40 : 41 : 42 . Although 
the discussion of our current publication is aimed towards 
quantum systems, our formalism may be applied with no 
change to many of these classical systems to rigorously 
derive, in the continuum limit, the well known classical 
glide motion (see Fig.QJ) occurring in the absence of in- 
terstitials/vacancies and high order effects. This rigorous 
result concerning the linear regime of continuum limit 
elasticity complements older works addressing detailed 
climb diffusion in various systems, e.g4£. We further re- 
port on new generalized glide equations in the presence 
of both dislocations and disclinations. As we illustrate in 
detail here, the basic physical ingredient leading to these 
generalized glide equations in solids is mass conservation 
which strictly restrains the dynamics. 

Recent years saw the extension of investigations of 
defect dynamics in classical media to quantum sys- 
tems^ addressing electronic liquid crystal and other 
phases 44,45 in which the electronic constituents favor, in 
a certain parameter regime, the formation of an ideal 
crystal like stripe pattern which may then be perturbed, 
through a cascade of transitions, to produce a rich variety 
of phases. Such stripe patterns are observed in the high 
temperature superconductors and other oxides* 6 - and in 
quantum Hall systems 47 . Defects naturally alter the local 
electronic density of states allowing for spatial (and tem- 
poral) inhomogeneities of electronic properties 48 . Follow- 
ing the general notion that melting occurs by the conden- 
sation of topological defects, e.gm&, we may naturally an- 
ticipate that the study of topological defects is pertinent 
to the understanding of quantum phase transitions be- 
tween various zero-temperature states. Much unlike clas- 
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sical physics where statics and dynamics are decoupled, 
in quantum systems, time and space are deeply entan- 
gled. As a consequence, we expect the zero point kine- 
matics of defects to be of paramount importance both 
to the character of the zero-temperature states and to 
the nature of the phase transitions. Only quite recently, 
along with Mukhin, two of us 4 ^ presented a thorough 
analysis of the problem of topological quantum melting 
of a solid in 2+1 space-time dimensions. At first sight, 
this follows the pattern of the famous Nelson-Halperin- 
Y oun g49j50 t neor y f classical melting in 2D. However, 
the quantum problem is far richer. Leaning heavily on 
the formalism developed by Kleinert 9 , it was shown in 44 
that the transition from elastic solid to nematic (or 'hex- 
atic') liquid crystal is closely related to the vortex (or 
'Abelian Higgs') duality in 2+1D. The nematic quantum 
fluid can be viewed as a Bose-condensate of dislocations 
subjected to a 'dual' Higgs mechanism. In this formal- 
ism, the rigidity of the elastic medium is parameterized in 
terms of gauge fields ('stress photons'). In the quantum 
fluid, the shear components acquire a Higgs mass due to 
the presence of the dislocation superfluid. As it turns out, 
such a substance is at the same time a conventional su- 
perfluid, which may now be viewed as an elastic medium 
which has lost its capacity to sustain shear stress. 

The glide principle as known from metallurgy amounts 
to the observation that dislocations only move easily in 
the direction of their (vectorial) topological charge. In 
the nematically ordered state, the Burgers vectors of the 
dislocations forming the condensate are oriented along 
the macroscopic director leading to an anisotropic screen- 
ing of the shear stress^. Here, the dislocations coast 
freely only in the direction of their director. Accord- 
ingly, the elastic propagator can only be fully screened 
in the direction perpendicular to the director. Even 
more striking, it turns out^i that in the quantum field- 
theoretic setting the glide principle acquires a meaning 
which goes well beyond the conventional understanding 
of this phenomenon. In fact, the main goal of this paper 
will be to further generalize these notions beyond the lin- 
earized continuum limit o^, including a study of theory 
in higher dimensions and the incorporation of far richer 
defect configurations. 

In many standard texts, e.g&LS, the explanation of 
glide takes little effort. A dislocation corresponds with a 
row of particles (atoms) 'coming to an end' in the mid- 
dle of a solid. One way to move this entity is to cut the 
neighboring row at the 'altitude' of the dislocation and 
consequently move over one tail to cure the cut (Fig. 
The net effect is that the dislocation is displaced. This 
easy mode of motion is termed 'glide'. Moving in the 
orthogonal direction is not as easy. Let us try to move 
the dislocation 'upward'. This requires loose particles to 
make the row of particles longer or shorter ('interstitials' 
or 'vacancies') and since loose particles are energetically 
very costly while they move very slowly (by diffusion) 
this 'climb' motion is strongly hindered. Estimates on 
'climb' diffusion rates are provided in e.g~. Climb is 
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FIG. 1: An edge dislocation in a two dimensional medium. 
Here, allowed (matter conserving) cutting and reconnecting 
of 'atomic columns' gives rise to the glide motion of the dislo- 
cation line (at the center of figure) in the horizontal direction. 
By contrast, an emission or absorption of vacancies and inter- 
stitials is required for a climb motion (vertical motion of the 
dislocation in this figure) . Thus, an absence of additional ma- 
terial (mass conservation) strictly prohibits the climb within, 
continuum limit, linear elasticity. 



hindered to such an extent in real life situations (e.g. 
pieces of steel at room temperature) that it may be ig- 
nored altogether. The lower dimensional glide motion 
of dislocations is reminiscent of dynamics in the heavily 
studied sliding phases 51 in which an effective reduction 
of dimensionality occurs. 

Although this simple pictorial explanation suffices for 
practical matters in metallurgy, it is immediately obvi- 
ous that there is more to it. Dislocations are completely 
specified by their topological charge, the Burgers vec- 
tor. Expressed in this language, glide is a prescription of 
the dynamics of defects as dictated by their own topol- 
ogy: dislocations can only move in the direction of their 
Burgers vector! How can this be? 

As we review in section II, this kinematic role of topol- 
ogy arises most naturally in the context of field theo- 
retic description of plasticity. Linear elasticity does not 
capture interstitial/vacancy defects. The only degrees of 
freedom which are admitted by continuum limit linear 
elasticity are phonons (or, in dualized form, the stress 
photons). Given the absence of interstitials, the secret 
behind glide is that normal crystals are 'non-relativistic'- 
solids cannot displace along the time direction. As the 
field theory lacks knowledge of the interstitials, this con- 
dition translates (as we will illustrate in section IV) di- 
rectly into the glide constraint for the dislocation matter! 

From the glide constraint it immediately follows that 
gliding dislocations do not occupy volume. Volume is 
exclusively associated with interstitial matter. Accord- 
ingly, dislocations do not communicate with compres- 
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sional stress and in the dislocation condensate pressure 
stays massless because of the glide constraint! It better 
be so, because the dislocation condensate corresponds 
with a conventional superfluid, and conventional non- 
relativistic fluids carry sound. Amusingly, glide needs 
a preferred time direction and in a relativistic medium 
there is no such thing as a preferred direction. As Klein- 
ert and one of us pointed out^, since glide cannot be 
defined in a truly relativistic medium, dislocations have 
to couple to compressional stress with the effect that in 
the relativistic quantum nematic crystal sound also ac- 
quires a mass: this incompressible state turns out to be 
nothing else than the space time of 2+1D general rela- 
tivity. 

The main limitation oS^iS is that they focus is ex- 
clusively on the special case of dislocation glide in 2+1 
dimensions within the realm of the linearized theory. Our 
main result, presented in section IV, amounts to expres- 
sions for new glide-type constraints in arbitrary dimen- 
sions. We further investigate what these constraints im- 
ply for both dislocations and disclinations via an anal- 
ogous study of the more fundamental double curl de- 
fect densities. In addition, we show in section V how 
these wisdoms come to an end in the full non-linear the- 
ory where dislocations can in principle climb, acting as 
sources and sinks of interstitials. 

This paper is organized as follows: After reviewing 
the basics of elasticity and introducing some notations in 
section II, we derive expressions for the dynamical cur- 
rents for an arbitrarily dimensional medium in section 
III. In section IV we deduce a general relation describ- 
ing the glide constraint as it acts on dislocation, discli- 
nation, and defect currents according to the linear the- 
ory, which is applicable in arbitrary dimensions. This 
allows us to analyze the particulars of dislocation glide 
in 3+1D (as well as in arbitrary higher dimensions). In 
high dimensions, dislocations are no longer particles but 
instead extended higher dimensional p-branes with p=l 
in 3+1D and p=2 in 4+1D. We will show that the glide 
constraint influences the center of mass motion of the 
branes while the relative (transversal) motions of these 
extended manifolds remain unconstrained. This amounts 
to a precise mathematical description of the 'restricted 
climb' known from metallurgy. We subsequently turn to 
the disclination and defect currents and presents proofs 
for deep- but also disappointing traits of these currents: 
the mass conservation underlying dislocation glide turns 
into simple conservation laws. Another important prop- 
erty of the constraint is represented in the symmetry of 
the constrained current. By identifying this constraint 
as applying only to the unique (rotational symmetry in- 
variant) 'singlet' component of the dynamical dislocation 
current, the relation of the glide constraint to compres- 
sion is firmly established. In the final section (Section 
(V)), we lift this to the full non- linear level incorporating 
the presence of a lattice cut-off and interstitial matter in 
the field theory formalism. 



II. ELASTIC ACTION AND DEFECT 
DENSITIES. 

Deformations in ideal crystals are parameterized by 
displacement fields u = R — R , where R and Ro re- 
fer to the actual position and the position in the ideal 
crystal of the constituent, respectively. The action of an 
elastic medium contains strain, kinetic, and external po- 
tential components. These are, in general, functionals 
of the displacements u and their derivatives. The ki- 
netic energy density term is simply K, = ^p(dtu) 2 . The 
incorporation of the kinetic energy density K, in the to- 
tal Lagrangian density leads to the variational equations 
of motion for the dynamics of u in both quantum sys- 
tems (with bosonic displacements u) and simple classical 
ones. Much throughout, we will perform a Wick rota- 
tion (wherein the real time t replaced by —ir) to Eu- 
clidean space-time. In the resulting imaginary-time Eu- 
clidean action, the Lagrangian density now becomes an 
"energy" density given by the sum of kinetic and poten- 
tial parts. The elastic portion of the Lagrangian den- 
sity follows from the gradient expansion in terms of the 
displacements. To leading order ('first order elasticity'), 
the elastic energy density is diU a Cij a bdjU b with Cijab the 
elastic tensor whose general form is dictated by the sym- 
metries of the medium. In the remainder of the article, 
we will employ Latin indices to refer to spatial compo- 
nents, while Greek indices will be reserved for space-time 
components. Einstein summation convention will further 
be assumed for repeated indices. In full generality, the 
action will further contain higher order derivatives and 
anharmonic strain couplings: 

So = J drdx £ [u, du,..., d rn u], (1) 
corresponding to the partition function 

Z = J Vu e -So[u,9u,..., d m u]^ ( 2 ) 

In the absence of singularities, the action above cor- 
responds with the simple problem of acoustic phonons. 
However, when defects are present, the theory becomes 
far richer. In its full non-linear form, elasticity is closely 
related to Einstein gravity^. In the 1980's, Kleiner^ 
achieved a considerable progress by his recognition of 
the underlying gauge field-theoretic structures, employ- 
ing the mathematical machinery of gauge theory to pen- 
etrate deeper into this subject than ever before; much of 
this is found in his book. This work focused primarily on 
the 'plasticity' of 3D classical media, and only very re- 
cently it was extended to the problem of quantum melting 
in 2+1D 44 . 

The key notion introduced by Kleinert is the dualiza- 
tion of the action Eq.JTJ into stress variables to recog- 
nize subsequently that stress can be expressed in terms 
of stress gauge fields or 'stress photons'. Topological de- 
fects then take the form of sources for the stress photons. 
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FIG. 2: Weingarten theorem: encircling a defect line in a 
three dimensional medium results in a displacement that de- 
pends exclusively on the starting/ending point x±. The dis- 
continuity occurs at the Volterra cut (rastered plane). 

This follows the same pattern as the vortex (or Abelian 
Higgs) duality where the phase modes of the superfluid 
are dualized in U(l) gauge fields, mediating the inter- 
actions between vortex sources. This analogy is quite 
close when only dislocations are in the game. The quan- 
tum nematics addressed by two of us earlier— are dual 
superconductors in the same sense that the quantum dis- 
ordered superfluids (Mott-insulators) can be regarded as 
dual superconductors, the difference being that in the 
dislocation condensate shear acquires a Higgs mass while 
in the vortex condensate super-currents acquire the Higgs 
mass. When disclinations come into play, we need to em- 
ploy rank two tensorial gauge fields ('double curl gauge 
fields'). On this level, the correspondence with grav- 
ity become manifest. As this construct is not widely 
known, we summarize the basic dualization steps in Ap- 
pendix A. 

There is a famous theorem due to WeingartenS which 
states that the singular displacement recorded by encir- 
cling a defect line in a three-dimensional elastic medium 
may always be expressed as a sum of a constant vector 
and an antisymmetric operator acting on the radial vec- 
tor x (see Fig. [2} 

j> (dl • V)u = u(x+) - u(x_) = b + &f(x - x ). (3) 

The reference point Xo can be arbitrary, but when fixed 
on the defect line it defines the Burgers vector b or 'dis- 
location charge' of the defect. The three-dimensional an- 
tisymmetric operator can be represented as the cross 
product of the radial vector x with the pseudo-vector 
f2: stfyi = O x x. The pseudo- vector f2 is the Frank 
vector or the 'disclination charge' of the defect. This 
theorem makes immediately clear that dislocation and 
disclination charges are related. For instance, a dislo- 
cation can be viewed as a bound state of a disclina- 
tion and an anti-disclination, while the disclination cor- 
responds with a 'stack' of an infinite number of disloca- 
tions. When addressing linear elasticity, the proper topo- 
logical current turns out to be a particular combination 
of dislocation- and disclination currents (the 'topological 
defect current—). 



Employing the dualization recipe of Appendix A in 
three classical dimensions, we can straightforwardly de- 
rive differential expressions for the dislocation and discli- 
nation densities respectively, 

af = e ijk djd k u a , (4) 
0" = ^e ljk e abc djdkdbu c , (5) 

where e denotes the 3D antisymmetric (Levi-Civita) ten- 
sor. For a defect line in a given spatial direction i, the 
charges in Eqs.|@J[5J|, derived in Appendix A via a dual- 
ization procedure, coincide with the Burgers and Frank 
vectors, respectively. To see this, consider, for instance, 
the edge dislocation of Fig. 0. Here, in going around the 
dislocation point there is a net change of the displace- 
ment fields. Couched in standard mathematical terms, 
the Burgers vector signaling the total variation in the 
displacement field is, as well known, 

b = be x = <j> Vu x ■ dl = J J V x Vu x ■ dA. (6) 

In Eq.JU), e x denotes a unit vector along the horizontal 
direction and dA is a planar area element within a region 
T> containing the dislocation point. We recognize the ar- 
gument of the last integrand in Eq. as the dislocation 
density of Eq.(@J- V x Vu x = a x ej. Similar considera- 
tions lead to the identification of Eq.JJJ as the angular 
mismatch (the disclination) density. A nice construct for 
its visualization is the Volterra cuA Needless to say, the 
derivation of the defect densities via integrals such as that 
in Eq.© is topological; these forms are valid for any crys- 
tal regardless of its fundamental constituents (whether 
they are classical or quantum matter of one particular 
statistics or another). In what will briefly follow, within 
the arena of dynamical defects (Section III), we will ele- 
vate the densities of Eqs.JUEJ m t° the zero (temporal) 
component of Euclidean space-time defect currents. As 
a curiosity, an old famous anecdote concerning time like 
defects is provided in— 

As stated earlier, dislocations and disclinations are not 
independent entities: a dislocation may be viewed as a 
disclination- antidisclination pair while a disclination is 
an infinite stack of dislocation lines. Thus, the densities 
of Eqs.JU 01 harbor redundant information. This re- 
dundancy may be removed by fusing the dislocation and 
disclination densities into a more fundamental (double 
curl) topological defect density, 

Vi = \£i 3 k£abcd 3 d b (d k u c + d c u k ) = 

Within this publication, most of our studies of these 
densities will be within the confines of the linearized the- 
ory of elasticity. When higher derivatives of the displace- 
ment field u appear in the action Eq.JTJ, the defect densi- 
ties discussed above will fail to capture all of the relevant 
singularities. The role of higher order corrections will be 
briefly touched upon in section V. 



Depending on the alignment of the charge vector and 
the three dimensional defect line, dislocations are clas- 
sified as edge (perpendicular, i ^ a) or screw (parallel, 
i — a), while the disclinations may be of the wedge (par- 
allel, i = a), sway (perpendicular, i ^ a) or twist vari- 
ants. 



III. DYNAMICAL DISLOCATION- AND 
DISCLINATION CURRENTS IN ARBITRARY 
DIMENSIONS. 

In this section, we will generalize the gauge theoretical 
description of defect currents to higher dimensions limit- 
ing ourself to the linearized level (first order elasticity) . 

As a defect drifts through the medium, the defect 
charge transforms into dynamical defect current. For the 
particular case of the 2+1-dimensional medium, such is- 
sues were addressed for the first time by two of us^. As 
a first extension, let us now consider how these currents 
look like in 3+1D. One has now to dualize the full ac- 
tion including the kinetic term. This is a straightforward 
extension of the 2+ ID case, 
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e^u\pd\d p u a , 



with the local rotation given by 



. a /3 _ 1 



(8) 
(9) 

(10) 



The currents are now two forms in the 'lower' space time 
labels, actually keeping track of the fact that the defect 
lines spread out in world sheets (or 2-Branes) in 3+ ID 
space-time. Obviously, by taking one of these labels to 
correspond with the time direction one immediately re- 
covers the static densities, = a- and T™ = Of. The 
currents with two lower spatial indices correspond with 
the dynamical currents: the lower indices ij refer to a 
current in the i direction, of a defect line extending into 
the j direction. 

In solids (and other condensed matter systems), the 
medium cannot displace in the temporal direction r: 
u T = 0. Consequently, in the definition of Eq.©, 
the'upper' Burgers-labels a of the dislocation currents J° 
are purely spatial. Formally, these spatial directions can 
be regarded as 'flavor' degrees of freedom^. We derive 
Eq. © by employing a generalization the static current of 
Eq.© along with an application of the Weingarten the- 
orem (Eq.©). The local rotation operator Ea. (|10|) used 
in the definition has also been generalized to deal with 
the in-equivalence of space and time. When one of the 
indices a, ft is temporal, Ea. (|l(J|l reduces to the usual 9 
spatial rotation lj ti = aA The corresponding disclina- 
tion current then carries the Frank charge of the defect. 
What is the physical meaning of the operator Ea. Hl(J|) 
when both indices a, f3 are spatial? Due to the constraint 
on temporal displacements (u T = 0), the only label that 
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FIG. 3: Discontinuity in the velocity field or the 'velocity 
dislocation': such a discontinuity, though physically possible, 
violates the assumptions of the Weingarten theorem Eq.lJ^J. 
Depending on the orientation of the Volterra cut with respect 
to the 'velocity Burgers vector' it represents slip of surfaces 
(parallel, full line) or adding/removing material (perpendic- 
ular, dashed line). For simplicity, we use the word 'slip' for 
both motions. 



can take the temporal value t is p and the generalized 
rotation represents the velocity field 



,ab 



2 £abcV 



(11) 



The corresponding dynamical disclination current T*~ 
records discontinuity, a 'dislocation', in the velocity field. 
Physically, this discontinuity represents slip of two sur- 
faces relative to each other, acting as a 'defect factory' 
in the crystal (Fig. ©. However, the Weingarten the- 
orem, which is the underneath the definitions Eq. © - 
©, assumes that all the symmetrized strains u p ^ = 
i(9 p M <T + d a u p ) are smooth everywhere. Choosing one 
of the indices p or a to correspond with time and the 
other to space, this strain becomes a local rotation, 



(12) 



A nonzero value of the disclination current would 
mean discontinuity in the velocity field, which in turn 
violates the Weingarten theorem. To avoid confusion be- 
tween real disclination densities and these surface slips 
we redefine them as follows, 



rpa 

pv 



2 £ [ivXk^ abc 

d\d K d b u c 



(13) 



so that only the Frank vector is represented by the upper 
label (disclination 'flavor'). When generalized to 3+1D, 
the dualization procedure summarized in appendix A in- 
dicates that the surface slips which we left out will even- 
tually play no role, because they do not couple to stress 
photons. 
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Let us now derive the form of the topological currents 
in arbitrary dimensions. The key lies, once again, in the 
Weingarten theorem of Eq.(PJ which is valid in any di- 
mension. In D dimensional space, the Burgers vector is 
a _D-dimensional while the disclination charge is a ten- 
sor of rank D-2. The rank of the disclination charge 
may be ascertained from the fact that the antisymmetric 
tensor srf can be written as a contraction of Levi-Civita 
symbol and another, D-2 dimensional antisymmetric ten- 
sor: = £yfe 1 fe 2 ...fc D _ 2 felfe2 ... feD _ 2 . Thus, two, three 
and four dimensional disclinations are characterized by 
a Frank scalar, vector and antisymmetric rank 2 tensor 
respectively. The D+l-dimcnsional extension of the def- 
initions of the currents Eq. I|8I13|) are 

^i-to-i = E ^...^D-iuxd u d x u a , (14) 

-^Ul...jt4D-l = 2 £ A I i -/ J D-lf \£a 1 ...a D -- 2 bcdvd\dt,U .(15) 

Both of the currents of Eas. (|14ll5(l are antisymmetric 
in the lower indices, and represent oriented p = D — 1- 
branes. The disclination currents are antisymmetric also 
in their upper indices where the Frank charge is dis- 
played. We reiterate that within linear elasticity the 
true (double curl) topological defect density of Ea. l(TB|) 
removes the redundant information given by the inter- 
related dislocation and disclination currents. Extending 
the duality procedure of linearized elasticity (appendix 
A) to arbitrary dimensions, we easily deduce that the 
fundamental double curl defect current within a D+l di- 
mensional medium is 

V\V2...V D -\ _ 

'lr) 1 r) 2 ...r)n-i 

v D - lP *d K d p {d x u a + cU K ). (16) 

This current is antisymmetric under exchange of any 
two lower or any two upper indices. This defect current is 
related to the dislocation and disclination currents Eq. (fT4l 
- IH>j) via 

D-l 

V1..MO-! _ \ ^ (\i+lx rpV 1 ...l/ i - 1 V i+1 ...l>n-l i 

f /m...np-i ~ / A > "™,-i )11 ..., lD _ 1 t 

i=l 

'^"2 £ ^—f'D-iKx9 K J l/1 _ , VD _ 1 + 2(D — 2)! £ ^i---A 1 -D-itA 

y -£vi...VD-ip& e p\t,l---£D-2 ^K"^Ci...£d_2 ■ (17) 

To avoid surface slips, one of the upper indices in the 
definition Eq. (jl6|l must be temporal when addressing the 
disclination currents T . When all the upper indices are 
spatial, such components record only the derivatives of 
the dislocation currents of Eq. l(Ti|) ■ 

IV. THE ELECTRICALLY CHARGED MEDIUM 
AND GLIDE. 

Let us now switch gear, to derive mathematical ex- 
pressions of the glide constraint in terms of the currents 
introduced above. We consider an electrically charged 



medium ('Bosonic Wigner crystal'), both because it is 
interesting on its own right 55 56,57 , and also because it 
provides us with a convenient vehicle for the derivation of 
general expressions for the glide constraint. Later on, we 
will independently derive these glide constraints by direct 
mass conservation without resorting to local gauge in- 
variance to implement it. As mass conservation pertains 
to a scalar quantity, such a conservation law translates 
into a condition on a linear combination of the topological 
defect currents which is necessarily invariant under spa- 
tial rotations. Conservation laws (equivalent to gauge in- 
variance) may greatly restrict the system dynamics lead- 
ing to an effective reduction in the dimensionality. We 
will now illustrate how this indeed transpires in solids: in 
linear elasticity, mass ('charge') conservation allows only 
a glide motion of a dislocation. 

We proceed with the treatment of a charged uniformly 
charged medium by employing the standard electromag- 
netic (EM) gauge field formulation. As usual, in Eu- 
clidean space, the electrical currents defined in terms of 
the constituent particles are minimally coupled to the 
electromagnetic potentials via ij^A^. The spatial cur- 
rent components relate to the velocity of the charged 
particles, so that every particle with charge e contributes 
to the current as (e d-rU 1 ). For a medium with n e of such 
particles per unit volume, the corresponding current den- 
sity will be ji — (n e e) d T u l . The time components de- 
scribe the coupling between the Coulomb potential A T 
and the charge density j T . For sufficiently small strains, 
the density is [(n e e) (1 — diU 1 )], subtracting the constant 
contribution (compensating electrical background). On 
this linearized level, we should add the following term to 
the Lagrangian of quantum elasticity to describe the cou- 
pling of the EM field to the electrically charged medium, 

Cem = i(n e e) \Aid T u l - A T diu l \. (18) 

The electromagnetic fields further have their own dy- 
namics, described by the Maxwell term Luaxweii = 
jF^F^ with field strengths F M „ = d ll A v - d v A p . 

Ea. i|18fl must be invariant under EM gauge transfor- 
mations A^ — > Afj_ + d ll a(x v ) 1 with a(x v ) being an arbi- 
trary, non-singular scalar function. The Maxwell term is 
automatically gauge invariant but the demand of gauge- 
invariance of the minimal coupling term Eci. (|18fl implies 
that the currents are locally conserved. After partial in- 
tegration, 

5C E m = ijpdpu — > -iad^j^. (19) 

As a(x^) is arbitrary, it follows that the current is 
conserved, = - the standard result that gauge 

invariance implies the conservation of the gauge charge. 
Let us now see how this works out considering instead 
the displacement. Performing a gauge transformation on 
Eq. (|18|) one finds immediately a gauge non invariant part, 

SCem — i{n e e) [diad T u l — d T adiU l ] 

— > i(n e e)a(x' M ) [d T diU l - didru 1 ], (20) 
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and gauge invariance implies that 

[drdiU 1 - d l d T u l ] = 



(21) 



Although this derivation is very simple, this equation 
is a key result of this paper. Combining it with definition 
of the dislocation current in arbitrary dimensions Eq. 114(1 
and using the contraction identity, 



--Tafii . 



(D-l)\ 



0~tv S t \ 



(22) 



it follows that the requirement of conservation of electri- 
cal charge has transformed itself into a constraint acting 
on the dislocation current, 



ja 

■id-1 u ix...%r>- 



0. 



(23) 



This is none other than the glide constraint acting on 
the dislocation current in arbitrary dimensions! 

At first sight this might appear as magic but it is easy 
to see what is behind this derivation. In order to derive 
Ea. (|18fl from Eq. (|19fl one has to assume that the gra- 
dient expansion is well behaved, i.e. the displacements 
should be finite. This is not the case when interstitials 
are present, because an interstitial is by definition an 
object which can dwell away an infinite distance from 
its lattice position. Hence, in the starting point Eq. i|18|) 
it is implicitly assumed that the interstitial density is 
identical zero. The gauge argument then shows that 
gauge invariance exclusively communicates with the non- 
integrability of the displacement fields, Ea. (|21|l . These 
non-integrabilities are of course nothing else than the 
topological currents - the glide constraint is a constraint 
on the dislocation current. If the glide constraint was 
not satisfied, electrical conservation would be violated 
locally, i.e. electrical charges would (dis)appear spon- 
taneously, as if the dislocation is capable to create or 
destroy crystalline matter. In the absence of interstitials 
this is not possible and, henceforth, dislocations can only 
glide. The key is, of course, that, by default, dislocation 
currents are decoupled from compressional stress in the 
linear non-relativistic theory. 

We may, indeed, equivalently derive Eq.J33J) without 
explicitly invoking EM gauge invariance to arrive at 
Eq. (|21|l . Instead, we may directly rely from the very 
start on mass conservation- the continuity equation of 
the mass currents (which, as alluded to above, is equiva- 
lent to local gauge invariance) 



drP + V-i = 0. 



(24) 



To see how this is done directly, we compute the vari- 
ous mass current components j^. By simple geometrical 
considerations, within the linear elastic regime, the mass 
density 



P = Po[l - diu 1 }, 



(25) 



with po the uniform background value: the divergence of 
u (signaling the local volume increase) yields the nega- 
tive net mass ('charge') density variation at any point. 



Similarly, the spatial current density 
j = p d T u. 



(26) 



Compounding the mass continuity equation of Eq.JSJJ) 
with the physical identification of the current fEos. (|25l 
EU), we obtain Eq.J2U from which Eq.(|23J) follows. In 
section (fV)l , we will return to such a physical interpreta- 
tion of the glide constraint from this perspective in or- 
der to determine corrections to the glide principle which 
follow from anharmonic terms. We emphasize that the 
mass conservation law leading to Eo. (|21|l trivially holds 
in any medium regardless of the underlying statistics^ of 
potential quantum systems or their dimension. Further- 
more, within the linear elastic regime, such a discussion 
highlights the validity of Eo. (|21|l and the ensuing glide 
equation of Eq. (|23[) (when interpreted as density matrix 
averages) in crystals at any temperature in which strict 
linear order mass conservation condition is imposed on 
all configurations. Needless to say, as temperature is el- 
evated, a departure occurs from such an imposed linear 
order condition through the enhanced appearance and 
diffusion of interstitials and vacancies leading to climb 
motions. The restriction on the dynamics in this regime 
is captured by a higher order variant of Eo. (|23[l (Eq.J5U|)) 
which will be derived later on. 

Let us now pause to consider what Eq.J^HJ) means 
physically. In two spatial dimensions, Eq. l|23|) implies 
that the dislocation currents have to be symmetric^*: 
£abJa = Jx ~ Jy = 0- Now, consider Fig. (1). Here, the 
Burgers vector is pointing in the horizontal x-direction 
implying that = 0, while the glide constraint reduces 
to Jy = 0. This current (Jy) is, by its very definition, 
the climb current perpendicular to the Burgers vector. 

In three and higher spatial dimensions, the story is less 
easy, the reason being that the constraint on the motion 
is less absolute. This is of course known in the classic 
theory^, but the reader might convince him/herself that 
making use of Eq. I|23(l the analysis is much helped as 
compared to the rather pain staking effort based on the 
'intuitive' arguments. In 3+1D, the constraint of Eq. i|23|) 
becomes 



^-rabcJbc ^* 



(27) 



Let us first consider a screw dislocation. These corre- 
spond with dislocation currents of the form J" (i.e. the 
static fi = t component corresponds with the orienta- 
tion of the dislocation loop being parallel to the Burg- 
ers vector). It follows immediately that the constraint 
Ea. (|27|) is not acting on screw dislocations and, hence- 
forth, screw dislocations can move freely in all directions. 
Edge dislocations are the other extreme, corresponding 
with dynamical currents of the form Jfi where a and b 
are orthogonal. The condition /i = a corresponds with 
glide: an edge dislocation with its loop oriented in a di- 
rection (b) perpendicular to the Burgers vector (a) can 
still move freely in the direction of the Burgers vector 
(a). The displacement of a dislocation along the line is 
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vector lies in this plane, this brane is analogous to a screw 
dislocation in the sense that its motions are not affected 
by the glide constraint. The other extreme is the 'edge 
dislocation brane' with a Burgers vector perpendicular 
to the brane, say in the x direction. In this case the non- 
trivial currents correspond with Jy Zfl - As noted earlier 
for three-dimensional defects, topological currents do not 
record any motion taking place within the defect-brane 
but only in directions perpendicular to it. The bottom 
line is that besides the static current (density) Jy ZT1 the 
topological dynamical currents are Jy ZX and Jy ZW , rep- 
resenting dislocation glide and climb respectively. The 
glide constraint Eq. (|23l) forbids the latter 'in integral 
form', allowing climb of a certain brane element only at 
the expense of the brane volume taken by a neighboring 
element. 



FIG. 4: Restricted climb in the 3+1D medium: although the 
dislocation line cannot displace as a whole in the direction 
perpendicular to both itself and its Burgers vector (climb di- 
rection), local segments may climb at the expense of their 
neighbors' height. 



not a topological object and the current with two identi- 
cal lower indices (p = b) vanishes. The glide constraint 
only strikes when all three labels are different. Let us 
consider a dislocation line extending in z direction with 
Burgers vector in the x direction (see Fig. 0J. The only 
nonzero components of the current are J z/1 and the glide 
constraint becomes J zy = 0. This automatically forbids 
any motion in the y direction, that is perpendicular both 
to the dislocation line and its Burgers vector. The con- 
straint is 'leaky' due to the extended nature of the de- 
fect. The material needed for the climb of one segment of 
the dislocation can be supplied by an adjacent segment. 
The glide constraint has therefore only a real meaning 
through its integral form, 







dVe Ta b c J; 



a b ' 



(28) 



As illustrated in Fig. only the dislocation's 'center 
of mass' is prohibited to move in the climb direction. 
Local segments of the line may still move at expense of 
their neighbors, effectively transporting matter along the 
defect line. The 'leaky', locally defined, constraint of 
Eq. (|28fl corresponds with the intuitive idea of 'restricted 
climb' found in the elasticity literature. 

Having a mathematical definition of the glide con- 
straint at our fingertips enables us to address its in- 
carnations in even higher dimensional systems where di- 
rect visualization is of little use. Notwithstanding that 
such crystals are, of course, not to be found in standard 
condensed matter, higher dimensional glide constraints 
might have implications for 'emergence theories' of fun- 
damental phenomena resting on elasticity theory«^S£. 
Let us for example look how the constraint Ea. (|23|) acts 
in a 4+ ID crystal. A dislocation is now a 2-brane, say, a 
plane extending in y and z directions. When its Burgers 



The glide constraint and motions of general 
defects. 



Up to this point we focused on dislocations, the tra- 
ditional scenery for glide motions. However, as alluded 
to in section II, dislocation and disclination currents are 
not independent. The truly fundamental topological cur- 
rents are the (double curl) defect currents of Eq. . A 
natural question to ask is what the form of the glide 
constraint is on these fundamental entities. This task is 
made easy by Ea. l21|) . In what follows, we will derive 
new generalized glide equations for a medium in which 
both dislocations and disclinations are present. 

Let us first consider the most trivial case, the 2+ ID 
crystal. Excluding surface slips, Ea. l|17|l becomes 



SurTp — \d\ [e\^ a Jy + £\vaJ^ ~ £ pva,J\\ ■ (29) 



In two spatial dimensions, we immediately note from 
Ea. (|29() that the glide constraint implies that rj\ = 0. 

The extension to arbitrary dimensions is straightfor- 
ward. We start with the definition Eq. i|16|) and then em- 
ploy the Levi-Civita contraction identity to obtain 



^/^ai . . .ao- ipA^rai . . .an- 1 be 1) 



[JLT 3fib 3fic 
3pr & pb & pc 
$\t $\b $\c 



(30) 



The terms containing u T vanish. The no-slip condition 
(d fl d u v a = d u dfj,v a ) needs to be imposed supplanted by 
the glide constraint in its strain form of Ea. (|20|) . The 
outcome is that the scalar spherical tensor component of 
the general defect current tensor of Eq. I|17|l vanishes 



ai02---OD-l _ A 

'Iaia2...an-l 



(31) 



We next investigate what the glide constraint implies 
for disclinations. The constraint of Eq. l|31|l does not shed 
much light on this issue. Linear elasticity directly ad- 
dresses dislocations: as we already argued, one of the 
upper indices in the defect current of Ea. (|16f) needs to 
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be temporal in order to record disclination currents. On 
the other hand, it is clear that disclinations somehow do 
know about glide because they can be viewed as an in- 
finite stack of dislocations. As an example, consider a 
wedge disclination: the material added by the Volterra 
cut is proportional to the Frank charge and this should 
not change over time. It turns out that the matter con- 
servation associated with a general defect distribution is 
related to the Volterra cutting procedure. After the cut 
has been applied, no additional material should be in- 
troduced, which is the same as the requirement that all 
symmetrized strains and derivatives thereof are smooth 
everywhere, including the locus of the Volterra cut. This 
condition is hard wired into the proof of the generalized 
Weingarten theorem Eq.J2J. As we will now show, the 
ramification of this principle for disclinations is a conser- 
vation law. 

To make headway, it is convenient to first consider Eu- 
clidean Lorentz-invariant D + 1 space time. The non- 
relativistic case will turn out to be a special case which 
directly follows from imposing the condition of the ab- 
sence of time like displacements (u T — 0). What follows 
rests heavily on elastic analogs of identities in differential 
geometry which are discussed in detail in the last part of 
Kleinert's boolA First, we introduce the tensor 

R^pa = - d u d p )dpU a (32) 

representing the Riemann-Christoffcl curvature tensor in 
the geometrical formulation of the theory of elasticity 
(the it's are crystal displacements). The smoothness as- 
sumptions underlying the Weingarten theorem may be 
expressed as, 

= (d^du-dud^dpu" + d a uP), (33) 
= {d^-dud^dxidpv? + d CT u»), (34) 

which in turn imply smoothness of the displacement sec- 
ond derivative 

= {d^du - d u dp)d x d p u a . (35) 

Cycling through the indices /i, v, and A of the smoothness 
equation Ea. (|35|l . we obtain the Bianchi identity for the 
displacement fields, 

dpRvX,po + d v Rx p , P a + d\Rnu,pcr — 0. (36) 

Starting from the other end, let us analyze a candi- 
date for a disclination conservation law, corresponding 
with d^Tfiviv^.'.'.v^LZ 1 ■ This can be expressed in terms of 
the Riemann tensor where a refers to a string of indices 
labeled by a, 

d ^ T n{v} = \ £ ^{u}K\£{ a }padpd K d\d p u a 

= ^£p{u} K x£{a}pad IJ ,(d K dx - dxd K + dxd K )d p u a 

— 2 ^ ^{v}K,X £ {a}padp 1 R K X,pa ~ ^^ p {v} ' 

The first term is zero as a consequence of the contrac- 
tion of the Bianchi identity Ea. (|36|l and the Levi-Civita 



symbol Ep{„} K x- This implies that relativistic disclina- 
tion currents are conserved. The non-relativistic case is 
just a special case: the vanishing of time like displace- 
ments means that all upper labels are space-like (a — > a) 
and it follows, 

In the 2+ ID medium, only wedge dislocations exist 
and the message of the conservation law Eq. I)38|) is clear: 
the Frank charge- the angle of an inserted wedge in the 
Volterra construction- behaves as a trivially conserved 
scalar component of a tensor. In higher dimensions, the 
disclination current similarly behaves as a conserved ten- 
sorial current. The defect density has information regard- 
ing both dislocations and disclinations. The disclination 
conservation law has separate ramifications for the de- 
fect density. The fundamental condition Ea. l|35l) can be 
directly rewritten into a conservation law for the defect 
density of a similar form as for the disclinations, 

Wr. 2 .^ 1 - 0. (39) 

This is not surprising as defect currents are propor- 
tional to to disclination currents. 

This completes the picture: the conservation of matter 
mandates that the 'proper' disclinations currents are also 
conserved. However, the 'handicapped' dislocation cur- 
rents are not conserved a-priori (disclinations form their 
sources) but they have to pay the price that they can 
only glide. 

B. Symmetry properties of the constraint 

The static dislocations and disclinations of higher di- 
mensional media are geometrically complex entities such 
as lines, sheets or d-branes with Burgers vectors and 
Frank tensors attached. When in motion, these branes 
sweep the additional time dimension. For instance, de- 
fects in two space dimensions are point-like particles 
turning into world-lines in space-time, in three space di- 
mensions they form loops spreading out in strings etc. 
Nevertheless, regardless of the embedding dimensionality, 
all of these 'branes' share an universal property: there is 
a unique direction perpendicular to the brane. This di- 
rection can be related to the dynamical defect currents by 
contracting it with the D — 1 dimensional antisymmetric 
tensor having one index set equal to time, 

1 7a 

(-D-1)! b rib 1 b 2 ...bn-i J bib 2 ...bo-l ~ 

d l d T u a - d T d lU a = U?, (40) 

isolating the perpendicular direction i. Instead of the 
dislocation current J, we could have used here also the 
disclination current T, but our interest in this subsec- 
tion will be in the former. The comparison of the tensor 
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Eq.^Ol with the glide constraint Eo. l(2Tj|) illustrates that 
the glide condition turns into a constraint on the trace 
U- =0. For rank 2 tensors, the trace is the only invari- 
ant tensorial component of the D-dimensional orthogonal 
group (O(D)), corresponding with all rotations and the 
inversion in the B-dimensional medium. It follows that 
the constrained current is the only 'singlet' (scalar) under 
the point group symmetries of the crystal (point group 
symmetries constitute a subgroup of 0{D)). The con- 
jugate degree of freedom must have the same symmetry 
and this can only be compression, the only physical en- 
tity being a singlet under 0(D). We have identified the 
fundamental reason that glide implies the decoupling of 
dislocations and compressional stress^. 

Apart from rotations, the glide constraint is also in- 
variant under Galilean space-time translations. Needless 
to say, the constraint does not obey Euclidean Lorentz 
(space-time) invariance as the time direction has a spe- 
cial status both in Eq. f2T)jl and Eq. lf2ll|l . The origin of this 
lies, of course, in the definition of the crystalline displace- 
ment u and its role in the minimal coupling Ea. (|18[) . The 
displacements are defined under the assumption that ev- 
ery crystalline site has an equilibrium position which im- 
plicitly wires in that their world lines extend exclusively 
in the temporal direction. The crystal defect currents 
(Eq.(P20 and the disclination conservation law of Eq.JSHJ 
are invariant under Galilean transformations while they 
do not respect invariance under Lorentz boosts. 



V. THE GLIDE CONSTRAINT, THE LATTICE 
CUT-OFF AND ANHARMONICITY 

Contrary to our rigorous ('glide only') result concern- 
ing the linear regime of continuum elasticity, in real crys- 
tals dislocations do climb (albeit at small rates). As dis- 
cussed earlier by two of us^, interstitial matter will be 
exchanged when dislocations collide and this process re- 
leases climb motions. In what follows, we will rederive 
the glide constraint, yet now do so within a fully general 
framework which will enable us to address the implica- 
tions of both (higher order) non-linear elasticity and the 
presence of a lattice cut-off. Higher order corrections 
to linear elasticity modify the original glide constraints 
of Eo. l|23|) giving rise dislocation climb. We will leave 
the detailed analysis of the non-linear problem to a later 
publication^. In what follows, we outline how the full 
non-linear theory (including finite lattice size effects) cap- 
tures higher order effects such as climb. 

To achieve this aim, we return to the mass continuity 
equation invoked earlier (in unison with Eqs.lJlllESl) yet 
now, by examining contributions of higher order deriva- 
tives of the displacement field, we exercise far greater care 
in examining its ramifications. The continuity equation 
of Eq.J2U implicitly assumes that the density and current 
fields are functionals of local Eulerian (distorted lattice) 
coordinates. On the other hand, the displacement, stress, 
and other elastic fields are functionals of substantial co- 



ordinates (i.e. the coordinates defined relative to the 
undistorted lattice coordinates)- the Lagrangian coordi- 
nate frame. As was briefly done earlier fEos. 1)251 126(1 ). 
we express the local density and currents of Ea. (|24|) in 
terms of volume and velocities as p = jr and j = 
with po the mass of the ideal uniform medium in a unit 
volume in the undistorted original medium (i.e. the ideal 
background mass density) . Following a distortion, a unit 
volume element of the original medium now occupies a 
region of volume V. With these relations in tow, Eq. l24|) 
reads 

d T V+{v- V)V = D T V = V(V-v). (41) 

This equation can be interpreted as a law governing the 
change of volume of the elastic medium: the change in 
volume (the derivative on the left hand side) is exclusively 
dictated by the motion of the boundaries- the gradient on 
the right hand side corresponds with a surface integral of 
the velocity field. This is just a reformulation of the same 
basic constraint: the conservation of mass (or electrical 
charge). Throughout this paper, mass conservation was 
the primary ingredient leading to the glide constraint. 
To invoke the mass continuity equation in the form of 
Ea. (|41|) . let us express the actual atomic coordinates BP 
in terms of the Eulerian coordinates (henceforth denoted 
by r l ) by employing the identity, 

(¥r)r = + E ^—r d >T\ (42) 

or 1 ( to! j 

m—l 

where a is a (lattice) cut-off scale. Simplifying, we find 
that Eq. (|41|l may be recast as 



d T V = V (— T ) T did T u>. (43) 
or 1 

This is an exact expression (entailing corrections to all 
order in gradients of u) detailing the glide constraint. 
Retaining the leading order contributions and employ- 
ing V = 1 + diu\ we recover the 'familiar' linearized 
glide constraint of Ea. H20|) . The exact glide constraint 
Ea. (|43() can also be derived using the coupling to elec- 
tromagnetism as generating functional as used in section 
IV. 

How may we generalize this analysis so that it deals 
with the non-linearities up to all orders? A key is pro- 
vided by Eq. (|18|) . Due to changes in the volume of the 
medium arbitrarily large lattice deformations translate 
into current densities j^. In addition, we should use 
substantial- instead of local coordinates such that the 
spatial arguments of the EM fields and strains no longer 
coincide. A crystalline constituent that was originally at 
position R feels the EM field at point R + u. Hence, the 
full non-linear generalization of Eq. 1|18|) is 

C EM = £j^[0 T u*(R)i4<(R + u) + j4 t (R + u)]. (44) 
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Using the procedure of Section IV, the exact glide con- 
straint is straightforwardly derived from this expression. 

What can we learn from this exact form of the 
glide constraint? Let us specialize to the simple 2+ ID 
medium. Assuming a vanishing lattice constant o, the 
volume of an elementary cell is simply given as V — 
det(Sij + diiii). Inserting this in the exact glide expres- 
sion Ea. il43|) we find that it simplifies, 



Q = e ab J b i {5 ia + d i u a ). 



(45) 



The displacement derivative diu a includes both regular 
and singular components. When diU a is small compared 
to unity, the linearized glide constraint of Eq.J^Hl is re- 
covered. 

Let us now analyze the conditions required for a dislo- 
cation in climbing motion to satisfy the generalized con- 
straint of Eq. (|45|l . As before, we set the Burgers vector 
of a dislocation to b = be x implying = 0. A climbing 
dislocation then moves in the y direction resulting in the 
dynamical current (of, is the position of the dislocation) 



Jy = b d T y 6[a 



(46) 



Assuming that the dislocation does not glide simulta- 
neously, the current Eq. I|46fl is the only nonzero dynami- 
cal dislocation current. The exact glide constraint is then 
satisfied, if and only if 



= 1 



d y u y 



(47) 



Clearly, so long as we keep the displacements small 
compared than the lattice size (effectively prohibiting in- 
terstitials and vacancies), the condition of Eo.l|47[l is im- 
possible to satisfy and climb motion is strictly forbidden. 
However, as soon as particles are allowed to leave their 
initial positions in the crystal (interstitial density devel- 
ops), d y u y becomes finite and climb motion is liberated! 

More precisely, the right hand side of the condition 
Eq. (|47|l represents the derivative of Euler to Lagrangian 
coordinates §^ and if vanishing, it means that two neigh- 
boring sites of the original crystal are, upon deforma- 
tion, occupying the same place. This is none other than 
an interstitial 'event'. We can thus interpret the 2+1- 
dimensional exact glide constraint of Eq. 1|45[) as formally 
proving the old metallurgical maxim that "climb mo- 
tion is permitted if and only if interstitials/ vacancies are 
present" . This climb will have the effect that the inter- 
stitial is absorbed by the dislocation. 

We reach the conclusion that higher orders in the dis- 
placements (non-linearities, anharmonicities) have the ef- 
fect of restoring the physics of interstitials. Another 
natural question to ask regarding the continuum limit 
is whether there may be a hidden dependence on the lat- 
tice cut-off a. To investigate this issue, let us see what 
happens to Eq. l|43|) when we only keep terms which are 
harmonic in the displacements, while imposing no condi- 
tions on the number of derivatives. These have to do with 



the volume of the unit cell. Let us specialize to a sim- 
ple hypercubic lattice with sides defined by the vectors 
(i = l,...,D), 

00 m-l 

u j (Ri + m) - u j (Ri) = 6 l3 + —T d ? u3 - ( 48 ) 

m— 1 

The corresponding volume is a determinant of the ma- 
trix Eq. 148|) and in the harmonic approximation only the 
diagonal elements remain, 



D oo 



(49) 



i^j—l m=l 



Using this expression for the volume, the exact equation 
of constraint (Eq. I|43[) ) reads 



f ■ ■ T a 

c Tan...iD-i 1 'i 1 ...i D _ 



adi 



(50) 



The remarkable fact is that it is possible to collect all 
higher order derivatives in a simple exponent. 

Eq. H50fl amounts to a double check on the non-singular 
nature of the continuum limit. It corresponds with a 
resummation to all orders and if there would be a non- 
perturbative surprise it should show up in this equation. 
One sees immediately that this is not the case: in the 
long wavelength limit (q — > 0), the spatial derivatives 
on the right hand side vanish recovering the linearized 
equation of constraint fEa.Q23[lL 



VI. CONCLUSION 

What has happened in this paper? This pursuit was 
driven by a bewilderment with the glide phenomenon. 
Surely, it is very easy to explain using the intuitive ar- 
guments found in the elasticity textbook (Fig.l). How- 
ever, we got puzzled when we discovered that it plays a 
rather central and outstanding role in the field theoretic 
formulation**. It is just not something that one adds by 
hand, but instead it is deeply wired into the structure of 
the theory, and it is responsible for the theory to have 
sensible outcomes like the 'prediction' that non relativis- 
tic fluids carry sound. 

In this earlier work 44 , it was assumed that glide was 
at work as a prerequisite to find sensible outcomes. In 
this paper we nailed down its origin: assuming that mass 
is conserved in the non-relativistic elastic medium, dislo- 
cations have to to glide as long as the 'interstitial' (non- 
topological) displacements are small compared to the lat- 
tice constant. In the linearized theory the glide constraint 
becomes exact. We learned in section V that this is to 
an extent an artifact of the field theory since there is 
room for interstitials in the fully non-linear theory. How- 
ever, one might want to view the field theory as a Pla- 
tonic form, a mathematical ideal, and it is just remark- 
able that topology, absence of symmetry (broken Lorentz 



12 



invariance) and non-topological conservation laws (mass 
conservation) conspire to produce a kinematic law deter- 
mined by a topological invariant. 

Our key results in Eas. (|2QI43|l translate mass conser- 
vation into the mathematical language of the topological 
defects. These are quite useful because they make possi- 
ble to address the workings of the constraint in circum- 
stances where the visualization 'method' of the elasticity 
books is of little use. We illustrated this by analyzing dis- 
location motions in arbitrary dimensions, demonstrating 
that this is about restricted climb ('mass transfer is free 
inside the brane'). This is surely of some practical value: 
we imagine that it might be useful in the complicated, 
real life situations encountered by metallurgists. It will 
be surely of importance in the generalization of the du- 
ality by Zaanen, Nussinov and Mukhir4^ from 2+ ID to 
3+1D, a task which is still to be accomplished and not 
unimportant given the fact that quantum nematics are 
more likely be found in a three dimensional world. A 
related issue is that it still remains to be demonstrated 
that the Lorentz-invariant version of such a 3+1D quan- 
tum nematic has to do with 4D Einstein gravity^. 

Why are disclinations not subjected to glide con- 
straints, while dislocations are? Using our formalism 
we could address this question in detail (section IV). In 
a way, the outcomes are surprising. Although we find 
that indeed the mass conservation does not impose kine- 
matic constraints on the disclinations, we do find that 
they imply a disclination conservation law. Although 
at first sight this might appear as obvious, it is actu- 
ally quite confusing. Usually topological conservation 
laws are rooted in the presence of an order parameter 
and are thereby associated with some generalized rigid- 
ity. In this way, dislocations are about shear rigidity 
while disclinations are associated with curvature rigidity, 
one of these beautiful facts becoming clear in the geomet- 
rical formulation. Volume and mass are associated with 
compressional rigidity. How can it be that this seem- 
ingly non-topological entity volume can still give rise to 
conservation laws pertaining to topological currents? 

We conclude with a speculation. As emphasized 
throughout, the basic physical ingredient leading to the 
derived constrained glide dynamics in solids was mass 
conservation. In a formal setting, similar albeit more 
restrictive volume conserving diffcomorphisms (parame- 
terized by time) may be directly examined via the al- 
gebra in two spatial dimensions and its extensions, e.gS. 
We suspect that there might be a more fundamental way 
of casting our relations by relying on the intricacies of 
such symmetries of space and time. 
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APPENDIX A: DUALIZATION AND DEFECT 
DENSITIES 

In this appendix, we sketch the dualization proce- 
dure for an action describing a three-dimensional elas- 
tic medium. We will then elaborate on the relation be- 
tween the dual fields and the defect densities. The ex- 
tension of our treatment to other dimensions (as well 
as to imaginary time) is straightforward. The disloca- 
tion/disclination densities of Eas. l|4l5|l harbor a redun- 
dancy which may be removed by considering a single 
symmetric tensor encapsulating information about both 
dislocations and disclinations- the double curl topological 
defect density ij of Eq. lltill . 

Linear elasticity is valid only when the displacements 
are minute as compared to the lattice size. Within this 
regime, the Lagrangian density of Eq.Q is 



— hdiU CijabdjU . 



(Al) 



This is a functional of displacements u a and their first 
derivatives dtu a . The elasticity tensor C must adhere to 
symmetries inherited from the underlying lattice. Thus, 
as dictated by the triclinical point symmetry group in a 
three dimensional medium, it may contain up to 21 inde- 
pendent constants. [In general dimension d, the elastic 
tensor may have up to ^d(d + l)(d 2 + d + 2) independent 
components.] The displacements u have generalized mo- 
menta conjugate to them- the stresses given by 



dC 



d(d lU a ] 



— Cijabdjl 



(A2) 



From the aforementioned symmetry of the elastic tensor 
C, an Ehrcnfcst constraint encapsulating a symmetry of 
the stress tensor (of = a l a ) rigorously follows. By a Leg- 
endre transformation, we may obtain the Hamiltonian as 
functional of displacements and stresses 

H[u a 7 of] = a*d t u a - C[u a 7 d iU a (u a , <)]. (A3) 

As we are after the Lagrangian in its dual form, we may 
first express the Lagrangian through the displacements 
and stresses exclusively via 



£dual = cr^diU a — Tt, 
leading to the partition function 



(A4) 



(A5) 



The dualization of the linear elastic action proceeds 
by a complete removal of the displacement degree of 
freedom. This is achieved by splitting the displace- 
ment field into the smooth and singular components 
u = u s .,„ + umv ■ This prescription generalizes the stan- 
dard Abelian Higgs duality, e.gS, to the elastic arena. 
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The singular displacement field umv harbors the defect 
density while the smooth field u sm is everywhere regular. 
Much as in other dualities^, the smooth part can be inte- 
grated out, leading a constraint on stresses (momentum 
conservation) 

a ( < - 0. (A6) 

By virtue of the Ehrenfest constraint, this conserva- 
tion law fEci. HAGI) ) is valid also when the contraction is 
performed with respect to the second index. This prop- 
erty enables us to express the stress field as a double curl 
of the symmetric two-form x (also known as the elastic 
gauge field£) 

of = £ijk£abcdjd b Xk,c- (A7) 

The physical field (stress) is invariant under the gauge 
transformation 

Sxk,c = dktc + d c &. (A8) 

We are nearly done. The Hamiltonian Ea. (|A3|) may 

now be expressed in terms of the elastic gauge fields 
X- The leftover term with the singular displacement 
afdiU a MV leads to (after integrating by parts) 

jCcoupi. = Xk,c [26£ + e abc d a a b k ] , (A9) 

with the final identification of the dislocation and discli- 
nation defect densities of Eas. (|4l5Jl . 



As usual (e.g., the famous minimal coupling (j^A^) 
in electrodynamics), such densities multiplying a gauge 
potential (here Xk,c) serve as source terms. In Ea. l|A9ll . 
the charges appearing in the elastic minimal coupling are 
topological. 

The redundancy of the defect densities needed in the 
coupling term Eq. (|A9|I is obvious: from nine indepen- 
dent topological defect densities (six dislocation density 
components af and six disclination density components 
0f , connected by three relations), only three linear com- 
binations couple to the physical degrees of freedom. 

To overcome the possible confusion, the coupling term 
may be treated differently so that the Ehrenfest con- 
straint is utilized at the very start 

<cU a = ±« + al)d lU a = <±(cU a + d a u l ). (A10) 

After the integration by parts, the elastic gauge field x 
becomes minimally coupled to the double curl topological 
defect density of Eq. @ 

Ccoupl. — Xk,cVk- (All) 

In spite of its usefulness in linear elasticity, the dou- 
ble curl topological density r\ is insufficient when higher 
order derivatives become important and new degrees of 
freedom are introduced. 
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